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Abstract 

We show that the XYZ spin chain along the special line of couplings Jy + JxJz + 
JyJz — possesses a hidden — (2, 2) supersymmetry. This lattice supersymmetry 
is non-local and changes the number of sites. It extends to the full transfer matrix of 
the corresponding eight-vertex model. In particular, it is shown how to derive the 
supercharges from Baxter's Bethe ansatz. This analysis leads to new conjectures 
concerning the ground state for chains of odd length. We also discuss a correspondence 
between the spectrum of this XYZ chain and that of a manifestly supersymmetric 
staggered fermion chain. 
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1 Introduction 

The solution of the zero-field eight-vertex model by Baxter is a landmark in the theory of 
exactly solvable systems. The seminal papers [HIU [3] [U [S] present a variety of algebraic 
and analytic tools to compute the partition function, the eigenvalues and eigenvectors of 
its transfer matrix. The model is still under active study: in particular, Fabricius and 
McCoy showed that at the so-called root-of-unity points the spectrum of the transfer 
matrix possesses degeneracies not easily explained by its standard integrability alone 
[5] [71 [5] [5] • It is natural to attribute them to the presence of extended symmetries, possibly 
elliptic generalisations of the sl2-loop-algebra symmetries discovered for the six-vertex 
model [inj[IIJ[I2]. 

The topic of this paper is the eight-vertex model at a particular root-of-unity point 
with an extended symmetry. We utilise the quantum XYZ spin chain, whose Hamiltonian 
commutes with the eight- vertex model transfer matrix |13| . More precisely, we study the 
Hamiltonian 
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with periodic boundary conditions along the special line of couplings 



JxJy + JxJz + JyJz = 0. (lb) 



While it was already noticed in Baxter's original works that along the line of couplings 
([T|) the ground state energy per site remains zero in the thermodynamic limit, its finite-size 
ground state was addressed much more recently. In [T3J [TS] Stroganov argued that ([TJ 
possesses exactly two zero-energy ground states for TV odd. When — Jy (so that 
A = Jz/ Jx = ^1/2), the resulting critical XXZ chain has been extensively studied over 
the past decade |inilI71IIHllISll2nilIIlllllll3]- The components of the ground state in 
the basis with Sz diagonal possess some remarkable properties. They display a variety 
of relations with combinatorial quantities, such as the enumeration of alternating sign 
matrices and plane partitions. Many of the early conjectures, such as sum rules for the 
components, were then proved with the help of techniques such as the qKZ equations, and 
combinatorial tools [Ml 123 123 HZ] . 

The fact that the ground state energy is exactly zero in the XXZ chain at A = —1/2 
for a finite odd number of sites was proved by exploiting a hidden supersymmetry 
[TOl I2U] . A ground-state energy of exactly zero is a common characteristic of theories 
with supersymmetry [25] • While quite a number of lattice models possess scaling limits 
described by field theories with supersymmetry |29j . only a few are known where the 
supersymmetry is explicitly present on the lattice. Here, an unusual feature is that the 
supersymmetry operator changes the number of sites by one. While unusual, it is not 
unheard of; similar operators were for example studied in a spin chain arising from the 
integrable structures in four-dimensional gauge theory |5U] . 

Remarkable properties of the zero-energy ground state persist along the entire line Q 
[211 132 , 33, 34,^351. Using the convenient parametrisation 

J. = 1 + C, - 1 - C, - (C - l)/2, (2) 

computer results indicate that the ground state at odd N can be expressed as polynomials 
in C with positive integer coefficients. Non-linear recursion relations were observed for 
several components of the zero-energy ground states of ([T]) [3TJ [321133] (see also [35]). 
These relations are described by the tau-function hierarchies of the Painleve VI equation 
[37] . Moreover, very simple expressions for one-point functions in the ground state such 
as the magnetisation were found by summing the expansion around the trivially solvable 
point C — >■ oo [35]. 

Since the ground-state energy remains zero along the entire line ([l]) , it is natural to 
expect that the supersymmetry persists off the critical point. This is true in the scaling 
limit, where the XYZ chain is described by the sine-Gordon field theory. This field theory 
possesses symmetries that can be related to the affine quantum group t/g(sl(2)) with zero 
centre |38j . When = —1, the quantum group contains the Af — (2, 2) supersymmetry 
algebra. This is precisely the value of q corresponding to the scaling limit along the line 
([T|). We thus refer to this line henceforth as the "supersymmetric" line. 

The purpose of this paper is to show that the supersymmetry remains exact on the 
lattice all along the supersymmetric line. The main focus in most (but not all [391) earlier 
studies was the ground state of the system, which has exactly zero energy for any odd 
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N. In this article, we study not only the ground states of this XYZ spin chain, but also 
consider its full spectrum. We show that for any value of in certain momentum sectors 
this model possesses a hidden exact symmetry which we call lattice supersymmetry, a 
lattice version of the well-known Af = 2 supersymmetry algebra [28] . These results provide 
a systematic generalisation of the XXZ results to the XYZ setting. 

We illustrate this symmetry by explicitly diagonalising the XYZ Hamiltonian along the 
supersymmetric line. Figure [T] shows the example of chains with A*" = 6 and N = 7 sites: 
in particular momentum sectors, the corresponding Hamiltonians sites have exact common 
non-zero eigenvalues. This is a well-known feature of supersymmetric theories: states with 
E > are doubly degenerate. Thus, in our case the common eigenvalues are candidates 
for supersymmetry doublets. This is a consequence of the existence of supercharges which 
change the number of sites N by one, and serve as intertwiners for the Hamiltonians of 
the corresponding chains. Moreover, two exact zero-energy ground states are found for 
odd A^, and they are potential candidates for so-called supersymmetry singlets (who do 
not have a superpartner). In fact, we will show that in the sectors with momentum zero 
for odd TV and momentum tt for even TV, the states organise themselves into quadruplets: 
an eigenvalue occurring first in the spectrum for TV — 1 sites, appears twice at N sites, 
and once at TV -I- 1 sites. The quadruplet structure hints at the existence of two copies of 
the supersymmetry algebra. Indeed by taking into account the symmetry under flipping 
all spins, we find a A/" = (2, 2) supersymmetry algebra on the lattice. 




Figure 1: Rescaled eigenvalues ej(C) = 4i?j(C)/(3 + C^) for the XYZ Hamiltonian with 
TV = 6 and TV = 7 sites along the supersymmetric line as a function of C- The plots are 
restricted to the subsectors with momentum k = n and respectively. The solid lines 
correspond to exact common eigenvalues for the two chains. Moreover, the plot for TV = 7 
shows the existence of exact zero energy ground states. 

The transfer matrix of the eight-vertex model commutes with the XYZ Hamiltonian, 
so the two have the same eigenstates. Thus one might expect that these properties carry 
over to the full model. As we will show, the existence of the Af — (2, 2) supersymmetry 
algebra is deeply rooted in the eight-vertex model, and can be derived from the Bethe 
ansatz. Moreover, the Bethe ansatz analysis leads to a novel characterisation of the ground 
states of the XYZ chain ([T| for odd TV. 

The XYZ chain is closely connected to another model with lattice supersymmetry. 
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introduced in 40 and thoroughly analysed in [Tnilll]|15]. It describes spinless fermions on 
a chain with nearest-neighbour as well as the usual on-site exclusion. Connections between 
the two at the critical point are already known; energy levels coincide with those of the 
XXZ chain at A = —1/2 with a momentum-dependent twist [19j . A similar correspondence 
also holds in the open chain with a magnetic field [20] . Moreover, the ground states of the 
fermion models display a variety of relations with the combinatorics describes above |43j . 
Therefore, it is natural to ask if a similar equivalence holds off the critical point. Indeed, 
several connections between the zero-energy ground states for chains of odd length and 
the ground states of a staggered version of the fermion models were observed in [3S1 HI] . 
In this work, we provide more evidence for the connection of the two models. We observe 
the systematic existence of common eigenvalues in their spectra after proper identification 
of the relevant parameters, even in sectors where the supersymmetry is not fully realised. 

The layout of this paper is the following. In section [2] we introduce supercharges for the 
XYZ chain which change the number of sites and study their properties. In particular, we 
show that Hamiltonian can be obtained as a quadratic form of the supercharges in special 
momentum sectors. The combination with various other symmetries of the XYZ chain 
will lead to a second lattice supersymmetry, and thus explain the quadruplet structure 
in the spectra. In section [3] we change our point of view and study the system using 
Bethe ansatz for the eight-vertex model. After recalling Baxter's original approach, we 
present a derivation of the supersymmetry from the Bethe equations, and show that it is 
a symmetry of the full transfer matrix along the supersymmetric line, connecting systems 
of different sizes. After this we proceed with a conjecture on the nature of the ground 
states for chains of odd length. The relation of the present model to the theory of fermions 
with nearest-neighbour exclusion is discussed in section [4] we observe that upon proper 
identification of the relevant parameters, the spectra of both theories have exact common 
eigenvalues. We suggest a mapping between the two models, based on the Bethe ansatz 
solution of the eight-vertex model. Finally, we present our conclusions and various open 
problems in section [5] Some technical details are relegated to appendices. 

Throughout the paper we report observations obtained from exact diagonalisation of 
small system sizes. Most of them show patterns which seem to hold for general N, and 
thus we formulate them as conjectures. 

2 Supercharges and the XYZ Hamiltonian 

In this section, we establish the relation between the XYZ chain (jlj and lattice super- 
symmetry. We start with a review of elementary symmetries of the XYZ chain and then 
proceed to the definition of the supercharges. The analysis of their properties under 
translation allows to show that they are nilpotent and write the Hamiltonian as quadratic 
forms of the supercharges. Next, we study the interplay with other symmetries such 
as parity and spin reversal, and establish the quadruplet structure of the eigenstates in 
special momentum sectors. 

In the sequel, we will frequently deal with operators that change the number of sites of 
the chain. Hence, we will indicate by a subscript the length of the chain the corresponding 
operator acts on. For example, we denote by Qjv the supercharge acting on the Hilbert 
space for chains with N sites etc. 
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2.1 Elementary symmetries 

We start with a bit of notation. We denote hy Hn — (C^)*^^ the usual Hilbert space for 
N spin-1/2 particles in a chain. We will use the standard orthonormal basis in which the 
operators (j| are diagonal {a'^ , a''^ , are the usual Pauli matrices). Its basis vectors are 
labeled by configurations a — aia2 ■ ■ ■ ctN with aj ~ + (spin up) or — (spin down) for 
the j-th spin, and the (t| operator acts according to 

(T||ai • • • Qfj • • • ajq) = OLj\a\ ■ ■ ■ aj ■ ■ ■ a^). (3) 

Let us indicate here some elementary symmetries of the Hamiltonian ([l]) , valid for 
any choice of Jx,Jy,Jz- First, it is invariant under translation. We will often use the 
translation operator Tjv acting on basis states oi "Hn according to the usual rule 

T^lai ■ ■ ■ a^-iaN) = |aArai • • • aN-i)- 

The translation invariance of the Hamiltonian implies that it commutes with this operator, 
and therefore they can be diagonalised simultaneously. Consider a state such that 
— tN\ip), then cyclicity implies (Tjv)^ — 1. Therefore the eigenvalue is an iV-th 
root of unity. Writing = e^'' we see that the momentum k has to be an integer multiple 
of 2tt/N. In this work, we focus mainly on momentum fc = for chains of odd length, and 
momentum k — n for chains of even length. 

Moreover, the Hamiltonian is invariant under reversal of the order of all spins. We 
thus define a parity operation P/y through 

-Piv|aiQ!2 • • • aN-ictN) = \o:nC(n-i ■ ■ ■ a2ai)- 

Obviously, we have = 1 and therefore eigenvalues ±1. 

Finally, notice that the Hamiltonian is a quadratic form of the Pauli matrices. Therefore 
it remains unchanged under global rotations by an angle n around any of the x—, y— or 
0— axis. Let us first consider the z axis. The rotation is given by Sn where 

N 

^^ = aiV|...a^ = exp|^(l-a|). (4) 

The right-hand side makes evident that Sn has the eigenvalue ±1 on configurations with 
an even/odd number of spins — . 

Considering instead rotations by the angle tt about the a;-axis leads to the conclusion 
that the Hamiltonian commutes with the spin reversal operator Rn defined through 

Rn = cr^a^-'-cr'^N- (5) 
The operators Sn and Rn have the following (anti-)commutation relation: 

SnRn — {~^)^ RnSn- (6) 

This implies in particular that for odd N the spin-reversal operator couples the sectors 
with even and odd number of spins down. Therefore any eigenvalue of Hn has even 
degeneracy and is at least doubly degenerate. 
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2.2 Supercharges 

We proceed with the definition of the supercharges. It is useful to start by recalling the 
usual N = 2 supersymmetry algebra [21] • It is built from two conjugate supercharges 
Q, that are nilpotent, i.e. their squares are zero: — {Q^Y = 0. A further symmetry 
generator is the fermion number F. The supercharges obey the relations Q] = Q 
and [F, Q^^] = — Q^. Hence Q increases the fermion number by one, while decreases 
it. The Hamiltonian is given as anticommutator H = {Q, Q^}. It conserves the fermion 
number and commutes with Q and . The fact that H is of this form implies that 
its eigenvalues are non-negative. Any states with zero energy are automatically ground 
states, and annihilated by both supercharges. Thus they each are a supersymmetry singlet. 
Conversely, all states with positive energy E > Q form doublets of the supersymmetry 
algebra, with the fermion numbers of the two states in a doublet differing by one. 

We here define the analogous supercharges on the lattice. We start with a generalisation 
of the work [2D] and introduce operators qj which map T-L^ to Hn+i through action on site 
j according to the following rules (the subscripts denote the positions of the corresponding 
spins) : 

qj\ai ■ ■ ■ Qfj-i + aj+i 
j 

qj\ai ■ ■ ■ - aj+i 

3 



We see that while states with spins + at site j are annihilated by g^, the spins — are 

transformed into pairs +- 1- or with weight 1 and —C, respectively. The pair creation 

implies a shift of the spin sequence ctj+i • • • by one site to the right. This includes 
a "string" (—1)-'^^ which is crucial in the following. To respect the periodic boundary 
conditions here, we need to build eigenstates of the translation operator. We thus introduce 
an operator that creates a pair of like spins ++ or on sites iV + 1 and 1, namely 

golai • ■ ■ aAr-i+> 0, (7c) 

N 

gol"! • • • oiN-i-) = - ( I + ai • • • UN-i + ) - CI - "1 • • • oiN-i - )]■ (7d) 

N \ 1 N+l 1 N+l I 

The operator go a-cts always on the last site irrespectively of length of the chain with no 
string attached. 

With the help of these "local" operators q^ we construct a supercharge Qat in the 
following way: suppose that |?/;) is an eigenstate of the translation operator Tjv with 
some eigenvalue tAr. Then we define Qn\'^) — unless % = (—1)^"^"^. If = (— 1)^+-^ 
however, we define 




For C = these operators decrease the number of — spins by one but map a state of 
definite total magnetisation to another state of definite total magnetisation [2D]. For finite 



aN) = 0, (7a) 
) = (-1)^"^ ( |ai • • • Oj-i + + ctj+i ■ ■ ■ aw) 

^ 3 J + 1 

- CI"! • • • "j-i - - "i+i • • • aw) ) • (7b) 

j .7 + 1 / 



• • • aN) = 
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C however, this is not the case because the two states on the right-hand side of Q differ 
in magnetisation by two. This is related to the fact that the XYZ Hamiltonian can flip 

pairs +- 1- and of adjacent spins, and thus conserves magnetisation only mod 2. 

Because of Qq the operator Qpf seems to distinguish the last site of the chain from the 
others. However, translation invariance removes this distinction, and we claim that Qn 
is a well-defined mapping between the momentum spaces with = (—1)^"*"^ . To show 
this it is useful to understand the transformation properties of qj under translation. From 
their definition it is not difficult to show that 

Tn+iQn = i~lfqo. (8) 

Consider now an eigenvector lip) of the translation operator for the chain with N sites: 
TnIiP) = tN\ip), tN = (—1)''^"'""^. Upon action with Q^q we produce a vector \cf)) ~ Qn\^) 
belonging to Hn+i- The application of T/v+i to this new vector leads to 



As tjv = (—1)^^^ the new vector is an eigenvector of TV+i, with eigenvalue —t^ — tjv+i, 
what proves our claim. Thus, we have T/v+iQ^rT^^ = —Qn- 

The fact that the operators Qat are bona fide mappings between the momentum spaces 
of interest is crucial to the supersymmetric structure of the XYZ chain Q that we describe 
now. First, the supercharges have "square zero" in the sense that 

Qn+iQn = 0. (9) 

Hence they can be thought of as fermionic. Second, if restricted to the subsectors with 
translation eigenvalue tjq = (—1)^+^ the XYZ-Hamiltonian can be constructed from the 
supercharges and their Hermitian conjugates. The latter are defined in the usual way: if 
is a vector in Hn, and \(p) in Hn+i then (?/;|Q]y|(/)) = (^IQatIV')*- With this definition, 
the Hamiltonian can be written as an "anticommutator" 

Hn = Qjv-iQjv-i + QnQn- (10) 

The proofs of ^ and ( [To| ) are elementary but cumbersome. We present the details in the 
appendices |A.l| and |A.2| Here, we study their consequences for the eigenvalue spectrum, 
and thus give an explanation of the common eigenvalues for systems of different size. We 
will see that all these properties are familiar from the theory oi J\f ~ 2 supersymmetric 
quantum mechanics |28j . 

First of all, consider the eigenvalue equation H^l^) — E\iIj). Projecting back on 
we find 

\\Q^j,_,m\' + \\QNH)r ^ E\\\m'- (n) 
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It follows that the spectrum is no n- negative: all £' > 0. Let us first concentrate on strictly 
positive energies -E > 0. For a chain with N sites these energies come in pairs in the sense 
that a given positive eigenvalue occurs in the spectrum at either iV + 1 or — 1 sites. 
This can be seen as follows. The structure of the Hamiltonian in ^ and (10) results in 
the commutation relation 

H-n+iQn — QnHn = 0. (12) 

Hence, if is an eigenvector of Hjy with eigenvalue E in the subspace with — (— 1)^+^, 
then QjvIV') is either zero or an eigenvector of Hp^+i with the same eigenvalue E. Likewise 
QIv-iI''/^) is either zero or an eigenvector of iJ^v-i with the same eigenvalue E. However, 
one of the two vectors Qn\^)^ <9]v-iIV') must vanislj^ Hence, every eigenstate with 
non-zero energy is part of a doublet 

m.QNm and so m,Q\^im- 



Conversely, zero-energy states are unpaired (singlets). From (111 it follows that they 
must be solution to the two equations 

QAr|^)=0, g]v_ilV^>=0. (13) 

It is well known that the solutions to these equations are related to the cohomology 
i^jv = kerQAr/ imQAr-i of the operator Qm- Every zero-energy eigenstate corresponds to 
a non-trivial (non-zero) element in i^^r. This can be seen indirectly: suppose that there 
are two linearly independent zero-energy states and |'02)i and assume that they are 
in the same cohomology class. This means that there is a state |(/)) for the chain with 
N — 1 sites such that 

|^l) = |V'2) + QiV-l|0). 

We act with Q\_i on both sides, and apply (13). This yields Q\_iQn-i\<I>) — 0, and 



by reprojection on |0) to ||(3Ar_i|0)|p = 0. This implies Qat.iI^) = and therefore 
IV'i) = l'02) ~ in contradiction to the assumption of linear independence. Hence every 
non-trivial cohomology class contains exactly one ground state. This has two consequences. 
First, every eigenstate with non-zero energy that is annihilated by Qjq can be written in 
the form ji/;) — QN-i\<i') where |0) is an eigenstate for the chain with A^ — 1 sites. We will 
frequently use this property in the next section. Second, the number of zero energy ground 
states is the number of distinct non-trivial elements in the cohomology. We formulate the 
following 

Conjecture 1. For odd N there are two non-trivial elements in S)m, whereas for even 
N there are none. 



iTo show this we write Hn = i?^' -|- i?^' with _ff^' = Qj^Qiv and H^^^ = Qiv-iQ]v_i- Then 
^N^^N^ = H^'h^^' = 0. Therefore, their respective eigenspaces associated with non-zero eigenvalues 
are orthogonal. If H^'li/i) = E\tp) for some E > then /f^'l?/)) = 0, by reprojection on \ip) we find 
IIQ]v-iI'/'>IP = 0' and therefore Q]v-iI'^) = ^- ^* follows that QnW is non-zero because otherwise \ip) 
would be a zero-energy state. Thus we found a doublet di/"), QjvIV"))- A similar argument applies to the 
eigenstates of and leads to pairs Q]v_iIV'))- Finally, as _ffjv commutes with and J^^' it 

follows that all non-zero eigenstates of our Hamiltonian organise in pairs in the sense stated above. 
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We checked this conjecture up to iV = 11 sites by evaluation of the row and column 
ranks of the rectangular matrices Qn- Notice that for N odd it implies Stroganov's 
conjecture |14j on the existence of two zero energy ground states, provided that one can 
prove that they occur in the zero- momentum sector. 

Before proceeding, let us make the following comment: both the nilpotency and the 
supersymmetric structure of the XYZ-Hamiltonian studied in this section are only vahd in 
certain momentum sectors. The nature of their derivation, given in the appendices | A . 1 1 and 
|A.2[ reveals that this restriction comes from matching the periodic boundary conditions. 
Everything else follows from local relations. Thus, we conclude that the supersymmetry 
has to be present in the full problem for chains of infinite length. 



2.3 Spin-reversal symmetry 

In this and the next section we provide a detailed discussion of the relation between 
supersymmetry and the other symmetries of the Hamiltonian introduced in section |2.1[ 
Let us give a motivation for this. The relation Qn+iQn ~ implies that we cannot relate 
chains with TV sites to TV + 2 sites by sole use of the supercharges defined in the previous 
sections. Yet, a detailed inspection of the spectra for small system sizes suggests such a 
connection. The most simple example is the the eigenvalue 3 + C^, appearing once in the 
spectrum for N = 2 sites, twice for N — 3 sites and once for = 4 sites. The inspection 
of other non-zero eigenvalues reveals similar patterns. This hints at a larger symmetry 
algebra which we analyse in this section. 

The existence of these degeneracies can be explained through a remarkably simple 
observation. The supercharges Qjv introduced in section |2.2| treat spins up and down 
in a very asymmetric way. However, the Hamiltonian commutes with the spin-reversal 
operator ([s]), as pointed out in section 2.1 Therefore it seems natural to introduce the 
spin-reversed version of Q^: 

Qn = Rn+iQnRn- 

In order to understand the implications of this operator on the spectrum we need to work 
out the algebra generated Qn, Qn and their adjoints. The full list of relations reads 

QnQn-1 = QnQn-1 = 0, Qn-iQn = Qn-iQn = 0' (14a) 

QjvQjv + Qjv-iQjv-i = 0, QjvQjv + Qjv-iQjv-i = 0, (I4b) 

QnQn-1 + QnQn-1 = 0, Qn-iQn + Qn-iQn = 0- (14c) 

The Hamiltonian can be written as anticommutator of either set of supercharges 

Hn = QnQn + Qn-iQm-i = QnQn + Qn-iQn^i- (14d) 

In fact, even any linear combination uQn + PQn with jap-f-l/Jp^lisan admissible 
supercharge that will generate the Hamiltonian. We omit the proofs of (14b I and ( 14c|) as 



Instead, we point out the striking analogies between (14) and the N = (2,2) super 



they are tedious and very similar to the proof of nilpotency outlined in appendix A.l[ 



symmetry algebra in two dimensional quantum field theory (see for example |45j . chapter 
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22). The latter consists of four supercharges Q±, Q±, Hamiltonian H, momentum P and a 
fermion number F. The algebra is defined through the relations 

= q1 = 0, (15a) 

{Q±,Qt}-0, (15b) 

{Q+,Q+} = A, {Q_,Q_} = A*, (15c) 

{Q+,Q_} = H + P, {Q+,Q_} = H-P, (15d) 

[F,Q±] = ±Q±, [F,Q±] = tQ±. (15e) 

together with the conjugation relations Qj_ = Qip and q1j_ = Qzp. The definition H and P 
implies that they commute with all supercharges and the fermion number. The operators 
A and A* are central elements. Non-zero values of the latter are usually an indication of 
topological sectors j46j which occur generically in non-compact spaces. 



The similarity between the two algebraic structures (14) and (151 is certainly not a 
coincidence. In fact, as mentioned in the introduction, the scaling limit of the lattice 
model is described by the sine-Gordon field theory at the supersymmetric point. Precisely, 
the limit is N ^ oo, C with g — (N^^"^ finite, and it yields the sine-Gordon theory 
with bare mass g and coupling /3 — -^/IGtt/S in conventional units. It was shown in 
that the sine-Gordon field theory possesses an affine quantum group symmetry [/g(sl(2)), 
q = exp(— Sivr^//?^) with zero centre. For /3 = -^/IGtt/S the latter is known to contain the 
Af — (2, 2) supersymmetry algebra [37], which indeed corresponds to our coupling. 

Let us establish a dictionary between the lattice and field theory quantities. The first 



three lines of (14 1 and (15 1 are in one-to-one correspondence, provided that we identify 
the field theory charges with some linear combinations of the lattice supercharges, and 
replace anticommutators by appropriate (and natural) expressions for the lattice operators. 
We see that there is no equivalent of the central charges A, A* in the lattice theory as the 
corresponding anticommutators are zero. We attribute this to the fact, that we have a 
finite, discrete and periodic space, and thus no topological/solitonic sector. The fourth 



line of ( 15 ) defines the field theory Hamiltonian and momentum. Whereas the former 
is related to the lattice XYZ Hamiltonian in a direct way, the identification between 
momenta in field and lattice theory needs to be supplemented by a shift by {N + l)7r 
mod 2tt. The need for a shift becomes even more plausible when taking into account that 
the lattice equivalent of fermion number appears to be the number of sites of the spin 
chain. This is consistent with the Jordan- Wigner transformation appearing in the analysis 
of continuum limit of the XYZ chain [351 . 

It is known that all non-zero energy states of a theory possessing the J\f ~ (2,2) su- 



persymmetry (15) are organised in quadruplets of the form (lip); Q+It/i), Q.jV'); Q-i-Q-IV")) 
with the state {ip) being annihilated by Q_ and Q_|_. The states in such a supermultiplet 
have all the same energy, and momentum, but differ in their fermion number. Given the 
similarities between the field theory and lattice algebra, it is natural to ask if there is a 



quadruplet structure in the lattice model. Indeed, we now show that the relations ( 14 1 
imply that it exists, at least in the momentum sectors considered in this paper. From 
the last section, we already know that the eigenstates {ip) of the lattice Hamiltonian 
Hn are part of doublets (IV')) Q]v-i IV')) or (IV')i Qa'IV'))- Without loss of generality, we 
focus on the second case (|V'), QtvIV'))- Let us consider the vector QnI^^)- There are two 
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possibilities: it may either be non-zero or zero. First, suppose that Qi^\'ip) is a non-zero 
vector. Clearly, it has the same energy as QatIV')) and we might wonder if they coincide. In 
fact, we show that this cannot be the case, and that they are rather linearly independent. 
If there is linear dependence, then there must be non-zero numbers A and ^ such that 
^Qn\iP) + P^Qn\iI>) = 0. We show that A = by applying Qjy from the left. For the 
first term, we use QjyQjvlV') = H]\[\ijj) ~ E\ip); for the second term, we make use of 
the anticommutation relation (14b) and write Q\^Qn\'4') — ^Qw-iQIv-i IV") = 0. Thus, 



we are left with — Ai?|-0) = but because of i? > we must have A = 0. Likewise, 
one shows that /i = 0. Thus, the two vectors are linearly independent. This implies in 
particular that in addition to Qjv_il'0) = we have the equation Qjv_ilV') = 0. Next, 
we increase once more the system size: consider the vector Qn+iQn\'>P) — ^Qn+iGnW- 
This vector is non-zero as again can be shown by an indirect proof: if for example 
Qtv+iQatIV') — then there would be some vector for the chain with N sites such 
that Qn\iP) = Qn\4>)- If we premultiply this relation by Qjy we find on the left-hand side 
QjvQjvlV') = ^Qjv-iQIv-iIV') = because of (14b). Thus, the right-hand side becomes 



Q^nQn\4>) — what implies Qn\4>) — 0. Yet, this is in contradiction to QnW 0, 
proving our claim. We cannot apply more supercharges in order to increase the length 
of the chain because the state Qn+iQn\iP) is annihilated by both Qn+2 (trivially), and 



Qn+2 (because of (14c)). Thus, we have constructed a quadruplet of one state at N sites, 
two states at TV -I- 1 sites, and one state -I- 2 sites, all of them having the same energy E 
with respect to the corresponding Hamiltonians: 



The preceding construction assumes that Qn]"^) is non-zero. Let us now consider the 
second case Qm]^^) =0. As > this can only be the case if there is a vector |(/)) for 
the chain with TV — 1 sites such that — Qn-i\4>)- Consider now the state QN-i\<t>)- It 
cannot be zero: otherwise, we could write = Qn(Qn-i\4>)) ~ —Qn{Qn-i\4>) = —Qn\'^) 
what contradicts our assumptions as we started from a doublet (IV')i Qa'IV'))- Moreover, 
the state Qis!-i\4>) is linearly independent from 1-0), as follows from the same argument as 
above. We have Qn\'^) = QnQn-i\4>) — —QnQn-i\4>) i a-nd hence a similar quadruplet 
structure as before, this time however with one state at iV — 1 sites, two states at N sites, 
and one state iV -I- 1 sites: 



[W.QN-M,QN-i\<t');QNQN-i\<t')) , Qn-M)^H)- (16) 

Therefore, all states with non-zero energy must be part of a quadruplet. We see that 
our argument leads automatically to a degeneracy in the "middle" of such a quadruplet. 
In fact, there is a non-trivial conserved charge that maps between these two states. It is 
given by 

Cn = Q\Qn = —Qn-iQ\^i- 

The anticommutation relations imply that it commutes with the Hamiltonian and has 
square zero: 

[Hm^Cn]^^, and = 0. 
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Moreover, its Hermitian conjugate is the "spin-reversed" operator Cj^ — RnCnRn- They 
have the character of fermionic ladder operators. Indeed, let us consider the quadruplet 



(16) containing two states = Qn-i\(I)) and \4>) = QN-i\(t>) at N sites. We find the 



following relations 



The other two states in the quadruplet, |</)) and QnQn-iI'P) ^ are annihilated by the 
corresponding operators Cn-i and Cn+i and their Hermitian conjugates. Thus, IV')) 



can be thought of a doublet inside the quadruplet ( 16 1 . 

Even though C^r is a bilinear in the supercharges, it still can be thought of as fermionic 
in the following sense. The symmetry operator Sjq defined in section [2. 1| anti-commutes 
with the fermion: CmSn + SnCn = 0, as can be shown using ([6| and the fact that 
Sn+iQn + QnSn — 0. If we now suppose that the state is an eigenstate of Sn with 
eigenvalue s = ±1 then the anticommutation relation tells us that \rp) is also an eigenstate 
of Sn, however with eigenvalue —s. Therefore, we see that the fermionic operators Cjv 
provide a mapping between the sectors with odd and even number of spins down. For 
chains of odd length, this conne ctio n is already established through the spin-reversal 
operator, as explained in section 2.1 Namely, the states {\tp), can be mapped onto 
each other through spin reversal, as the state \(j>) is an eigenstate of the spin-reversal 
operator Rn-i when A'^ is odd. For chains of even length, however, the spin-reversal 
operator fails to connect the two vectors, whereas the operator Cn does this independently 
of the number of sites. 



2.4 Parity symmetry 

In the last part of this section, we analyse the relation between the supercharges and 
the the parity operation. Using the definition of the "local" supercharges Q we find the 
simple transformation laws 

Pn+iQ] = {-l)'^^^qN-j+iPN; PN+iqo = QoPnTn- 

Let us now consider a state for a chain of N sites with both definite translational 
behaviour Tpflip) = tN\ip) and definite parity P/v|'0) — PNli^), Pn = ±1- The parity 
operation reverses momentum, as can be seen from the relation PnTj^P^ — T^^ . This 
implies that the translation eigenvalue tj^ must be solution to = 1. Obviously, this is 
compatible with = ( — 1)-'^+^. Applying our rules, we find 

Thus, we find that (1) for N odd, tjv = 1 the action of preserves parity and (2) for N 
even, = —I the action of reverses parity. 
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We studied the parity sectors by means of exact diagonalisation of the Hamiltonian up 
to iV = 11 sites. For odd N we observed that the spectrum of the parity odd sector is 
contained in the parity even sector. Thus, we arrive at the following 

Conjecture 2. For odd N and zero momentum, the spectrum in the odd parity sector 
Pat = —1 is contained in the spectrum of the even parity sector = +1. 

3 Supersymmetry and the transfer matrix of the eight- 
vertex model 

In this section, we consider the transfer matrix of the zero-field eight-vertex model. The 
main tool here is the Bethe ansatz for the transfer matrix established by Bax:ter [3J IH |S] . 
We provide a derivation of the supersymmetry from this point of view, generalising the 
result of |19j from the critical point to the entire supersymmetric line. 

After some basic definitions in section |3.1[ we provide a brief review of the Bethe 
ansatz for the eight- vertex model in the root-of-unity case in section |3.2[ in particular 
recalling the necessary change of basis of the Hilbert space. In section [3?3] we establish 
the supersymmetry in the new basis. To make contact with the supercharges defined in 
our previous discussion, we have to transform back to the canonical spin basis what is 
discussed in |3.4[ This leads to some new conjectures on the nature of the zero-energy 
ground states of the XYZ chain of odd length. 

3.1 Basic definitions 

We start by recalling elementary facts about the eight-vertex model on the square lattice 
135]. Each edge carries a classical Z2 "spin" variable ±, corresponding to occupied/empty 
or spin up/down. The configurations are restricted in such a way that each vertex has an 
even number of spins down: the eight allowed vertex configurations are shown in figure 
[2] We associate a Boltzmann weight to each vertex, and the weight of a given lattice 




aabbccdd 



Figure 2: The vertex configurations of the eight- vertex model (the bold edges have spin 
down). The associated weights a, 6, c, d are invariant under spin reversal in the zero field 
case. 

configuration is then simply the product over all the vertex weights. In the "zero-field" 
case, the weights are invariant under simultaneous reversal of all spins around a vertex. 
Thus, as shown in figure |2] there are four distinct weights, traditionally denoted by a, b, c, d. 
Suppose that the square lattice has say M rows and N columns, and is wrapped around 
a torus (periodic boundary conditions along the two directions). Then the model can 
conveniently be studied by the row-to-row transfer matrix Tn, whose matrix elements 
are defined as the sum over all configurations along a horizontal line, compatible with the 
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spin values on the vertical edges: 



{a'\TN\a) = E 
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The invariance of the vertex weights under spin reversal implies that [Tn-Rn] — 0. 
Moreover, the vertex rule implies that the transfer matrix conserves the number of down 
spins mod 2. Therefore we have [7~Ar, Sn] — 0. Conservation of the total number of down 
spins is only possible in the six- vertex limit = (or c = 0) . 

To proceed we parametrise of the vertex weights in terms of Jacobi theta functions, 
following the definitions of [HIl EH and [5U] : 

a = a{u) = p'&i(2r], q^)'di{u - rj, q^)di{u + rj, q^), 
b = b{u) = p^4{2rj, q^)di{u - 77, q^)di{u + r;, q^), 
c = c{u) = pi9i(2?7, q^)'di{u - ry, g^)i?4('" + 9^): 

Here w denotes the spectral parameter, 77 the so-called crossing parameter, and q the 
elliptic nome. Moreover, we choose the overall normalisation as p = 2/t?2(0, g)i?4(0, q^). 
This choice ensures that 

h{u)^a{u) + b{u) = di{u,q), (17) 

a function which we shall use quite often (the right-hand side follows from standard 
identities for Jacobi theta functions |50j). With this parametrisation two transfer matrices 
with different spectral parameters u, u' commute: 

[Tn{u),Tn{.u')]^Q. (18) 

This implies that the series expansion of the transfer matrix in the spectral parameter 
around any point yields a family of commuting operators. The most simple ones are the 
translation operator Tjv and the XYZ-Hamiltonian 

Tjv = /j(27/)-^rjv(r/), Hn - a{Ti)/b'{ri)TN{ur^T'i,{u)U^^. 

The Hamiltonian reduces exactly to our problem ([ij if the crossing parameter is set to 
rj = 7r/3. In this case, the variable C, used to parametrise the supersymmetric line is related 
to the elliptic nome through 

^ ^4(2-/3,9^); • ^''^ 

We are interested in using supersymmetry to study the eigenvalues and eigenvectors of 
the eight-vertex model transfer matrix. Because of ( 18 1 the eigenvectors do not depend 
on the spectral parameter, and so coincide with those of the XYZ Hamiltonian Hj^, up to 
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possible degeneracies. Such degeneracies do not seem to appear at generic values of the 
crossing parameter, but only at the special elliptic root of unity points r] = (mi7r+m27rT) /L, 
with mi, 1712, L integers and q = e'^"^, where additional symmetries are present [B]. In our 
case T] = 7r/3 we have already shown that in the momentum sectors with tjv = (—1)^^"'^ 
the eigenvectors organise into singlets, or quadruplets with the same value of E. Moreover, 
the eigenvectors in a given quadruplet each can be labeled by a distinct quantum number; 
two of them are for — 1 and + 1 sites, while we showed that the two at N sites have 
an even and odd numbers of spins down. Both the number of sites and the number mod 
2 of spins down are preserved by the eight- vertex model transfer matrix, so barring any 
accidental degeneracies, these correspond to distinct eigenvectors of the transfer matrix 
as well. Since Hn is obtained from the logarithmic derivative of Tn{u), the fact that 
their eigenvectors coincide makes it natural to hope that analogous structure occurs in 
the spectrum of the transfer matrix. We here show how at rj = -k the supersymmetries 
described above indeed extend to the transfer matrix, and so give relations among the 
eigenvalues. 

As an indication of the special properties occurring at rj = 7r/3, we note that the 
zero-energy states |\E'^) of Hm for odd A^, i.e. the supersymmetry singlets, have very 
simple transfer-matrix eigenvalues 7~Ar(M)|^'*) = Tn{u)\'^^) . They are given by [THITS] 



The simplicity of this expression stresses the special nature of the two eigenstates. In the 
sequel we will see that the study of the transfer matrix eigenvalues leads naturally to a 
distinction of these states from the other eigenvectors. 

3.2 Review of Baxter's Bethe ansatz 

Here we summarise the aspects of the coordinate-type Bethe ansatz [3J HJ [S] relevant to 
our derivation of the supersymmetry in the eight- vertex model at = 7r/3. 

Path basis. The transfer matrix of the eight-vertex model has no obvious particle-number 
conservation (such as conservation of the number of down-spins). This is a central difficulty 
when compared to the six-vertex model. In [4 Baxter developed a way to overcome this 
problem through the introduction of a basis upon which the transfer matrix acts in a way 
that resembles the six- vertex case. 

For A^ sites the new basis vectors are labeled by a sequence of integers £1,(2, ■ ■ ■ ,iN,^ n+i 
such that \ij+i £j\ — 1, j — 1, ■ . ■ , N . It is useful to think of a path starting at some 
height ii = £ with the restriction that consecutive heights differ by ±1. In the following, 
we will therefore frequently call the corresponding set of vectors in Hn the "path basis". 
A down step or particle occurs at site j if £j+i — £j = —1, and an up step occurs otherwise. 
The path is completely characterised by £ and the positions xi, . . . ,Xm of its m down 
steps. Hence for Xk < j < Xk+i the local heights are given by 



Twiu) = h{u) 



N 



^l{u,q) 



N 



(20) 
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(21) 



The basis vectors are given as an TV-fold tensor product of local vectors in 

N 

The factors are constructed from the local heights via 

|a>,+i,,) = 7?i(t+ (2^+1)77,(72)1+) + 7?4(f+ (2^+1)77,(72)!-), 

where |±) are the local spin— 1/2 basis vectors, and s and t arbitrary parameters such 
that the two vectors are linearly independent. 

The transfer matrix maps the set of these states onto themselves and conserves the 
number of particles, provided that the following condition is met for generic values of 
the crossing parameter rj the initial and final height are identical £1 = £n+i, and therefore 
N = 2m. In the special case of eUiptic roots of unity 77 = (toitt + ■m2TTT)/L however, this 
constraint can be relaxed to 



N — 2777 = Z/p, for some p € Z 



(22) 



because of the periodicity of the Jacobi theta functions involved in the construction of the 
vectors. The height difference between starting point and endpoint of the path is thus 
Lj) as shown in figure [3j Moreover, because of the periodicity of the theta functions it is 
sufficient to restrict the initial height to £1 = 0, 1, . . . , L — 1 in this case. 




Figure 3: A typical path for L = 3, two particles and y —\. 

It is instructive to compute the maximal dimension cIn oi the subspace spanned by 
these vectors in the root-of unity-case by simple counting of the paths. There are (^) 



arrangements of 777. particles, provided that (22 1 holds. Let us introduce an indicator 
function (^^(n) which is 1 if 77 = mod L, and otherwise. We have the convenient 
representation 



27rijn/ L 



L 

j=o 

We weight this by the number of arrangements and an additional factor L which takes 
into account the different choices for ^ = 0,1,. 1. Summation over m yields 

d,^L± SUN - 2777) . 2- g (cos ('-^y " 

m=0 ^ ' j=0 V \ / 
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For the case of interest L = 3,we find 

dN = 2^ + 2{-iy'^. 

We know that the dimension of the full Hilbcrt space Hn is 2^. If we assume that all 
the vectors associated to paths are linearly independent then we conclude that for even 
N the path basis is redundant. For odd however at least two vectors are missing, and 
thus the path basis does not span the entire Hilbert space. For small finite-size systems, it 
seems that exactly two vectors are missing, i.e. the existing 2^ — 2 vectors are linearly 
independent. We shall assume the linear independence in the following, and will conjecture 
later that the two-dimensional complement of the path basis at A'' odd is spanned by the 
ground states of the XYZ Hamiltonian. 



Eigenvectors and Bethe equations. The next step consists of decomposing the eigenvectors 
of the transfer matrix in terms of vectors in the path basis: 

L 

IV') = X]^^ ^■^{i\Xi,...,XmW,X\,...,Xm)- (23) 

e=i {xj} 

The summation over the positions of the particles is carried out in an ordered way 
1 < a^i < a;2 < • • • < Xm < N. Moreover uj is an L-th root of unity: = 1. The wave 
functions 'tp{t, Xi, . . . , Xm) are obtained through a Bethe-type ansatz. In order to describe 
it we need Baxter's "single particle" functions and "wave vectors" defined through 

g,(^,x) = e'^^- y"^+-^,-"-"^\ , e-^ = M^. (24) 
h{wi+x-2)h{'wt+x-i) h{uj - T]) 

Here we used the function 

we = is + t)/2-TT/2 + 2£rj (25) 

which is linear in £, and contains the free parameters. The numbers Ui, . . . , are the 
Bethe roots to be determined. With this notation the wave function is given in typical 
Bethe-ansatz form by 

tp{e\xi,...,Xm) = ^ A^5^(i)(f,a;i)fif^(2)(^-2,a;2)---S'^(„)(^-2(m-l),a;„). (26) 

Here the sum runs over all permutations tt of m objects. The m! coefficients satisfy 
the following relation: if r is a transposition exchanging j and j + 1, then we have for the 
permutation w' = n or the relation 

^ fe(M^(j+i) - u^U) + (27) 

An KUn{j) - UttO + I) - 2??) ■ 

The left-hand side is commonly interpreted as the (bare) scattering matrix between two 
particles with "rapidities" Un(j) and u^^j+i)- If they coincide then we find A^' = —A^, 
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implying that the Bethe wave function vanishe^ This will be very important in our 
analysis. 

The Bethe roots ui, . . . , Um remain to be determined. Baxter showed in [5j that if 
they solve the Bethe equations 



then (23l is an eigenvector of the transfer matrix T n{u)\^n) = Tn{u)\^Pn)- The corre- 



sponding eigenvalue can be obtained from the so-called TQ-equation 

TN{u)Qpf{u) = uj(I)n{u ~ r;)Qjv(M -I- 2ry) -|- + '>l)QNiu - 2ry), (29) 

where (f>N{u) = h{u)^ — t?i(u,q)^, and 

m 

i=i 

is an elliptic polynomial with zeroes at the Bethe roots. In particular, setting u = rj 
we find that the eigenvalue tj^ of the translation operator Tn is given in terms of the 
Q-function as 

- ^^'iTTT = n - ft (30) 



where in the last step we used (24 1. 



3.3 Derivation of the supersymmetry from the Bethe ansatz 

We now use the Bethe ansatz to establish for the case rj — Tr/3 the supersymmetry 
connecting systems with different numbers of sites N and ± 1. 



We start by noting that from (22 1, the number of particles m in the path when L = 3 
must obey 

TV - 2m = 3p 



^ Let us suppose that two Bethe roots have the same value, say Um-i = Um- Now we modify the sum 
over permutations in \2&\ according to = froir for some function / on the symmetric group 

Sm, where r is an arbitrary permutation of m objects. We choose t to be the transposition of m — 1 and 
m. For any tt define pre-images ni,n2 according to 7r(ni) = m — 1 and 7r(n2) = m, then we find 

i}{t\xi,...,x^)="^A^, Jl g^^^){l-2(j-l),Xj) 

X ff,r(n2)(^- - l).^ni)9,r(ni)(^ " 2(n2 - l),^:™^) 

According to our assumption we have A^/ = —At^. Moreover, from = and ( |24[ l, it is not difficult 

to see that gm-l((-, x) = gm.{i, x). Therefore we may write g^(„2) ~ 2(ni — 1), x„i )gn(ni)i^ ~ 2(^2 — 
l),x„2) = 9-n(ni){f'- 2{ni - l),a;„Jg^(„2){£- 2(n2 - l),a:,i2)- Using these facts, we see that 

. . . , Xm) = -ii{(\xi,. . . , Xm), 

and therefore the wave function vanishes. 
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for some integer p. This relation is compatible with the simultaneous replacement N — >■ 
N' ^ N — j , m ^ m' ~ m + j and p ^ p' = p — j for some integer j. The supersymmetry 
charge Qm studied in section [2] increases the number of sites by one. However, in the 
context of the Bethe ansatz it turns out particularly convenient to consider an action like 
that of Q\_i- we choose j = 1 and therefore decrease the length of the chain by one while 
adding a particle to the system. We discuss the relation with the supercharges studied 
previously in the next section. 

Our strategy is to construct from a given solution ui, . . . , Um of Bethe 's equations at 
N sites a new solution mi, . . . , Um+i at — 1 sites. We shall verify that a solution 
to this problem is simply given by Uj — Uj for j < m + 1 and Um+i — t^- Indeed, for the 
smaller system the first m Bethe equations with this choice become 

h{uj + f?) \ ^ ^ 2 TT ^{^j ^ + 2?^) h{uj — TT + 2ri) 

h{uj ^ 1]) J l^i'^j ^ "u-k ^ 277) h{uj — TT — 277) 



Using ( |28[ ) and the antiperiodicity h{u + tt) — --h{u), this equation reduces to 

h{uj — 77) h{uj — TT + 27;) 
h(uj + 7/) hiuj + TT — 277) 

It holds for generic Uj if 77 = tt — 277 mod tt, and thus in particular for the value r] = tt/S 
we are interested in. However, we still have to check the (tti + l)-th Bethe equation. We 
find 



On the right-hand side we recognise the eigenvalue of the translation operator for the system 



with N sites (30 1. We conclude that the operation is possible only if tjv = {—1)^^'^oj''^. 
But recall from the last section that for 77 = nm/L the number a; is an L-th root of unity, 
in our case thus uj^ — 1, and therefore we find the symmetry in the momentum sector 
with 

tAT = (-1)^+1. (31) 

Consistency thus requires a restriction to the momentum sectors studied in section [2] 

As a side comment, let us notice that we could have started from the TQ-equation 
with arbitrary uj, deduced the Bethe equations from the requirement that Tn{u) is an 
entire function and imposed the lattice supersymmetry. Asking for consistency would have 
led us to = 1. In the six- vertex limit where the working is completely analogous, write 
Lu — e"^ so that has the interpretation of a twist angle. We conclude that the twists 
leading to the symmetry here are (p = 0,±27r/3, as follows from the observations of [H]. 
These are precisely the values for which special ground state eigenvalues of the transfer 
matrix, as well as relations to alternating sign matrices, appear p!51 [T71 [T5] . 

As a next step, we determine the relation between the corresponding eigenvalues of 



the transfer matrix from ( 29 1 . For the Q- function we find 



Qn-i{u) = ]^ h{u - Uj) X h{u — tt) = —h{u)QM{u). 



k=l 
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Using this relation, we deduce that 



Tn{u) + h{u)TN-i{u) = 0. (32) 
Setting u — 7] we obtain a relation between the eigenvalues of the translation operators for 



both systems ^at^i = —t^ = (— 1) . This is consistent with (31 1. For odd N we obtain 



thus the zero- momentum sector (invariant under translation), whereas for even N it is the 
TT-momentum sector. This fits well the picture suggested by (24): the (to + l)-th particle 
with Um+i — has momentum km+i — tt, and therefore the eigenvalue of the translation 
operator is changed by a sign. 



Relation between eigenvectors. The preceding operation should manifest itself as an 
operation on the Hilbert space (or at least the special momentum sectors). In fact, we 
would like to introduce an operator Qn '■ Hn ^ 'H-N+i (not to be confused with Baxter's 
Q-matrix) such that the eigenvectors of the transfer matrix that can be obtained from 
the Bethe ansatz are related according to 

\lpN-l) = Ojv-ilV'Ar)- 

Twofold application of Qjy would lead to the injection of two particles with momentum 
TT. However, in this case the Bethe wave function vanishes. Hence we can write on the 
subspace spanned by the path basis 



0, or QnQn-1 = 0. 



We will now derive the explicit form of these operators, starting from the definition 



of the wave functions (26 1. To manipulate them, we need an explicit expression for the 



amplitudes A^^. In fact, their defining equation (27 1 can be solved up to a factor 



l<.i<j<.m 



2,/). 



(33) 



Let us consider the wave function (26 1 for to + 1 particles, one of them having momentum 
Um+i — It is a sum over permutations tt of {1, 2, . . . , to-|- 1}. For a start, let us consider 
in this sum only the permutations with 7r(r) = to + 1 for some fixed r = 1,...,to+1. Any 
such permutation can be decomposed according to tt = tt' o tt" where 7r'(TO + 1) to + 1 
and 



We have the signature sguTr" = (— l)™+i and thus sgnvr = (—1)'"+^ ''sgnTr'. We use 




this in order to evaluate the Bethe ampUtude (33) in terms of the permutation tt'. After 
some algebra one finds 



A^ = sgn7r' ]J Ku^.f^i) 

l<z<j<m 



m m+1 

^.'U)+2ri)\{h{u,-'n) n e'"'"'- 

i—1 j—r-\-l 
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We see that the only r-dependent term is the last product: a string of wave- vectors. Thus, 
we must understand how this affects the corresponding single-particle functions. For j > r 
we notice the identity 

e"'-''>9nU){i - 2i + 2, x^) = - 2(j - 1) + 2,x, + 1), 

which holds only because of ?7 = 7r/3. This is already enough to simplify the wave function. 
As tt' leaves m + I unchanged, we can think of it as a permutation of only m objects. This 
can of course be done for any value of r separately. Collecting the different contributions, 
we find after some algebra a recursion relation for the wave functions 

m m-t-1 



ilj{£;xi,...,x,n+i) =]^/i(wi~?7) X! 9m+i{£-2{r~ l),Xr) 



i—1 r—1 

X 1p{£; Xi,..., Xr-l,Xr+l + 1, ■ ■ . , X„i+i + 1). 

The wave functions on the right hand side are the ones involving only the Bethe roots 
ui, . . . ,Mm, thus precisely the ones for the problem at N sites. However, notice that the 
positions of the particles are not arbitrary: for the r-th term we have Xr+i -I- 1 — Xr-i > 2. 
The picture becomes a little more transparent if we consider the induced operation on the 
basis stated 

m+l Xr-2 

Q]v-il^;a;i,---,a;m)jv =C ^ ^ gra+i{£ - 2{r - 1) , x) 

r—1 x—x.r-i-\-l 

X \e;xi, . . . ,a;r_i,a:,a:r+i -^,---,Xm- l)7V-i (34) 
where we set xq = 0; C is a normalisation constant which can be chosen arbitrarily. We 



see that (5]v-i inserts a new particle at position x between existing ones at Xr-i and 
provided that they are at least two sites apart (in order to guarantee Xr-i < x < Xr — I)- 
More precisely the operation Q]v_i transforms locally two consecutive up-steps at (x, x-l- 1) 
to a single down-step at x, while shifting all particles on its right by one step to the left 
as illustrated in figure |4] 

This is weighted by Cgm+i{(i — 2k + 2,Xk = x). Notice that the local height at x -I- 1 is 
given by^a; = £ + x + 1 — 2{k — 1). Combining this with the definition of the single- particle 



wave function ( 24 1 , we conclude that apart from the string the weight can be expressed 
in terms of the local height between the two up-segments alone (to see this recall that 
h{u) is 27r-periodic and Wx+3 = Wx -\- 27r). More explicitly, we make the convenient choice 
C — 11^=1 ^(^f) s-iid find the local weight 

Cg™+i(£ -2k + 2,xk^ x) - {~\Yh{wt^^,f . (35) 



This completes the definition of the supercharges acting on the path basis. From (32 1, 
we conclude that they act as intertwiners for the transfer matrices: 

Tn{'u)Qn-i + h{u)QN-iTN-i{u) = 0, (36) 



•^Strictly speaking, QjY_i only acts on momentum states with tjv = (—1)'^"''"^. Clearly, these can be 
written as a superposition of path states. Thus, for simplicity, we present the action of the supercharge on 
the vectors in this decomposition. 
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Xk-1 X Xk Xk-1 X Xk - I 

Figure 4: Local action of Q\i_i- two steps up are transformed to a step down. 



on the momentum sectors with ijv — (—1)^^^. Of course, this relation only holds on the 
subspace spanned by the path basis. 

We finish this section by pointing out that - as in the case of the XYZ chain - the 
notion of "particle" is somewhat arbitrary. One could as well have chosen the steps up as 
particles. From the local vectors (21 1 we see that this corresponds essentially to exchanging 
the parameters s and t. This would lead to a second supersymmetry operation with the 
same local weights ( 35 1 , which transforms locally two consecutive steps down to a single 
step up, and thus resembling strongly the case studied in section |2] We will discuss their 
connections in the next section. 



3.4 Supercharges in the spin representation and the XYZ ground 
states 

Having found an operator Qjy (and thus Qn) defined through its action on states of 
the path basis it seems natural to ask how it acts on simple spin states, i.e. momentum 
states built from a spin configuration. Generically, the path states are rather complicated 
superpositions thereof, and hence we have to find the transformation relating the two 
bases. In order to work it out, we must address the question of incompleteness of the path 
basis for chains of odd length, pointed out in section [3?2l In fact, for odd N we have to 
define the action of Q}v-i missing two states which we denote by \^±). 

Let us first state a simple observation: as the \'^±) are not in the path basis, they 
cannot be obtained through the action of Qjy on any state in the Hilbert space Hn+i for 
the chain with A'^ + 1 sites. Second, we extend the definition of Q\^_i in the most natural 

way: QJv-iI^±) ~ ^- ^f course, the same reasoning applies to the operator Qn itself, 
and thus we have 

Qn\'^±) = 0, and l^-i) ^ Qn-i\^) for all \^) e Un-i- (37) 

Notice that this provides a consistent extension of the nilpotency property Qn+iQn — ^■ 
In a more mathematical parlance, we extend thus the definition of Qjq in such a way that 
the missing states are closed, but not exact with respect to the operators Qn- 

Given this extension, the next steps are to construct the relation between the path 
basis and the spin basis, and to understand the action of the operator Qjsi on simple spin 
states. We have not found a systematic construction for generic N . It seems that this is 
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related to quite non-trivial identities between Jacobi theta functions, as suggested by the 
most simple example Q2- We work out this special case in appendix |b] It suggests the 
following 

Conjecture 3. For any N the operator Q at is a linear combination of the supercharges 
Q N and Q n , defined in section [2| with coefficients depending on the free parameters s 
and t. The latter can be fine-tuned in order to make one of the coefficients vanish. As a 



generalisation of ( 36 1 we have the intertwining relation 

Tn{u)Qn-i + h{u)QN-iTN-i{u) = 

(and a similar relation for QN)on the momentum sectors with = (—1)^+^. 

It is easy to check that this equation is compatible with the commutation relations 



between the supercharges and the XYZ Hamiltonian (12). More importantly however, this 



conjecture ~ if true - has some interesting consequences for the ground states of the XYZ 



spin chain along the supersymmetric line. In fact, combining it with (37) we conclude 
that the two missing states |^'±) correspond to non-trivial elements in the cohomology of 
the supercharge Qn (or Qn), discussed previously. Therefore they are perfect candidates 
for the ground states of the XYZ Hamiltonian for chain with odd N. Indeed, we verified 
this conjecture up to A'^ = 9 sites by checking that the ground states obtained through 
exact diagonalisation of the Hamiltonian are indeed orthogonal to all states in the path 



basis provided that we impose the relation between C and the elliptic nome q given in ( 19 ). 
Thus we are led to 

Conjecture 4. For odd N the subspace of vectors that are orthogonal to the path 
basis is two-dimensional. It is spanned by the two zero-energy ground states of the 
XYZ-Hamiltonian ([T]) which are invariant under translation. 

This is perhaps the most surprising outcome of our analysis because it suggests that 
for r/ = 7r/3 the ground states at odd N cannot directly be obtained from Baxter's Bethe 
ansatz. Notice that this observation is different from the widely discussed question of the 
completeness of the Bethe ansatz (see e.g. [51 ) as here the Bethe ansatz (as it stands) 
does not apply to the missing states. 

These two states are thus eigenstates of the eight- vertex transfer matrix with eigenvalue 
A]y{u) = di{u, q)^ as was conjectured by Stroganov [HI US]. The conjecture was extended 
to the inhomogeneous eight- vertex model, defined by allowing on any site j a shift of the 
spectral parameter u ^ u — Uj. The conjectured eigenvalue is [34j 

Tn{u) = i)i{u - ui,q)'di{u -U2,q)- - ■ t9i(M - un, q). (38) 

The simple product structure of this eigenvalue suggests that there is a local mechanism 
leading to its existence. Here we seek to extend our conjecture 4 to the inhomogeneous 
setting. The construction of the path basis parallels the homogeneous case, with a slight 



modification of the local vectors (21 1. For site j they become 

^Ms+ {2£+l)v + uj,q^)\+) + Ms + (2^ + 1)V + wj, 9')|->, 

,0) 

-e+i 



,0) 

|ci>(i\ ^) ^Mt + m + 1)'7 - a') 1+) + ^?4(t + m + 1)77 - g2) 
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where |±) are the local spin— 1/2 basis vectors. Using this, we checked numerically for 
= 3,5 and 7 sites, and random choices for the spectral parameters Uj the following 

Conjecture 5. For odd N the subspace of vectors that are orthogonal to the inho- 
mogeneous path basis is two-dimensional. It is spanned by the two eigenstates of the 
inhomogeneous transfer matrix with eigenvalue ( 38 1 . 



4 Relation to lattice fermions with hard-core exclu- 
sion 

In this section, we present new observations about the connection between the XYZ 
model along the supersymmetric Une and the staggered supersymmetric fermion chains 
with nearest-neighbour exclusion considered in [33111] (see also |S2]). We provide a first 
step to construct a mapping between the models, based on the path description of the 
states in Baxter's Bethe ansatz for the eight-vertex model. Such a mapping was relatively 
straightforward to obtain for the XXZ case UHl [50] , because in both cases there is a 
conserved U{1) symmetry, and have closely related Bethe equations. Here there remains a 
U{1) symmetry in the fermion model, but there is no such manifest symmetry in the XYZ 
chain. Nevertheless, there appears to be a close relation between the spectra in the two 
cases. Moreover, we explain how there is evidence of the hard-fermion structure in the 
path basis of the eight- vertex model by exploiting the mod 3-periodicity of the heights in 
the path description. 



4.1 Conjectures relating the spectra 

Let us recall the model of |40II19I . describing spinless fermions on a periodic one-dimensional 
lattice with Nf^f) sites. The fermions are subject to the constraint that no two adjacent 
sites are both occupied. Defining the ordinary fermion creation and annihilation operators 
to be c| and Cj with {c^, Cj} = {cl, Cj} — and {ci, Cj} — Sij, the constraint amounts to 
restricting the usual fermionic Hilbert space to states annihilated by rijUj^i, where Uj = 
CjCj is the fermion number operator. Fermions respecting this constraint are annihilated 

and created by the operators dj = (1 — nj_i)cj(l — nj+i) and = (1 — nj^i)c^(l — nj+i). 
The model has explicit J\f — 2 supersymmetry: it is built from a supercharge 

^(/) 

Qif) = ^^4' 

where the Xj are non-zero real coupling constants (possible phases may be removed by 
simple gauge transformations. The Hamiltonian is given as anticommutator H(^fj = 

For periodic boundary conditions dj — dj^^^^^ on the fermions (known as 
Ramond boundary conditions), the Hamiltonian is 

^U) = H •^i^j+i(«'l+Wi + d]dj+i) + X]{1 - nj_i)(l - nj+i) . 
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This thus inckides a hopping term, a chemical potential, and a next-to-nearest neighbour 
repulsion. Notice that unlike in the XYZ chain where the magnetisation is not conserved, 
this Hamiltonian conserves fermion number for any values of the Xj. 

Following [331 SI] we now consider the case where length of the fermion chain is a 
multiple of three, and the coupling constants are staggered with period three. Then the 
problem is invariant under translation by three sites: if T(^f^ is the translation operator 
on the fermion chain then we have [H^f),T^j:^] = 0. In this case, one can show that for 
N(^f) = 3m the model has exactly two zero-energy ground states with m fermions in 
the "momentum sector" where T|?^j = (—1)"*"'""'^. The precise form of these ground states 
depends on the values of Ai, A2, A3. The most general case is analysed in [S^. Here we 
describe the choice Xi — y, X2 = I, X3 = y for some real y. This was the case studied in 
where we conjectured that after the change of variable 

C"^ = 1 + 8y2 (39) 

the two zero-energy ground states of the fermion chain at A'^(/) — 3to, and the two zero- 
energy ground states of the XYZ spin chain with N = 2m + 1 sites share some components 
which are polynomials in ^ and are related to a tau-function hierarchy associated with the 
Painleve VI equation [5T||5^|55]. 

It is natural to ask if the relation between the two models is deeper. Indeed, if we 
rewrite the spectrum of the fermion chain in terms of the variable C by using ( |39[ ) then a 
number of eigenvalues in the spectra of the XYZ Hamiltonian at = 2m and N = 2m + 1 
coincide exactly with eigenvalues of the fermion chain Hamiltonian ^zH^j-^ (the factor 4 is 
just an issue of normalisation). This statement can be sharpened by analysing different 
momentum sectors. As an example, we provide the characteristic polynomial det{E — Hn) 
for the XYZ Hamiltonian with 4 sites in the subsector with momentum k — n: 

{E - (C^ + 3)) {E - + 2C + 5)) {E - ie -2C + 5)) {E - 2 + l)) . (40) 



The characteristic polynomial det(e — ^(/)) for the fermion model at -/V(/) = 6 sites in the 
subsector with m = 2 particles and T^j^ = — 1 is given by: 

(6 - (1 + 4y2)) (e - (1 + 2y')) (e^ - (3 + 4y') e + 2{1 + 2y' + 2y')) . (41) 



If we set e — E/A and use the change of variables (39 1, then (41 1 coincides with (40 1 up 
to the factor E"^ and an unimportant global numerical factor. Hence we see that upon the 
change of variables, the spectra coincide with the exception that the zero-energy states 
are absent in the XYZ spectrum. This coincidence of the XYZ spectrum at iV = 2m 
on the sector with momentum vr, and the fermion model at A^(/) = 3m on the sectors 
with T^^-^ — (—1)™+^ appears to be systematic for small m, but different multiplicities of 
various eigenvalues occur for m > 4. Studying the spectra up to m = 6, we are led to the 
the following conjecture: 

Conjecture 6. The spectrum of the XYZ Hamiltonian i/jv for N = 2m sites in the 
sector with momentum tt coincides with the spectrum of the staggered fermion chain 
4:H(f) with N(^f) — 3m sites in the sector where T^j^ — (—1)™+^ if variables are changed 



according to (39), with two exceptions: (i) the eigenvalue E — is missing in the XYZ 
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spectrum and (ii) the two models lead to different multiplicities of the eigenvalues. 



This conjecture identifies sectors of the two models where their supersymmetries are 
exactly realised. Yet, it appears that the connection is even deeper. We analysed the 
relations between the models for antiperiodic or Neveu-Schwarz boundary conditions 
dj+Nf^fj = —dj on the fermions. In this case, the supersymmetry of the fermion model 
is broken. The spectrum needs no longer be positive, and indeed the ground state has 
negative energy. These boundary conditions are equivalent to a twist in the Hamiltonian, 
leading to the term —XN^f)^i{d\dN(f^ +'^jv(j:)'^i)- This sector is unlikely to share properties 
with the momentum sectors discussed so far in this paper because they have explicit 
unbroken supersymmetry. We found however coincidence with the spectrum of the XYZ 
chain of even length N = 2m and zero momentum. We illustrate it once again by showing 
the explicit characteristic polynomials for m = 2. The XYZ Hamiltonian for iV = 4 sites, 
and momentum k = has the characteristic polynomial 

(E-A) {E~iC- If) {E-iC + If) - {e + l)+2E + 3)' + 8 (C^ - l)') 

The characteristic polynomial of the fermion model at -/V(/) = 6 with Neveu-Schwarz 
boundary conditions restricted to the sector T|?^^ = 1 is given by 

(, _ 1) (,3 „ (2y2 + i-^+2e (2y2 + l)^ + 8y') {e' ~ e (V + l) + V) ■ 



Again, if we set e = E /A and perform the change of variables p9| we find that the two 
polynomials coincide up to some unimportant numerical factor. Studying small systems 
up to m = 6 we are led to the 

Conjecture 7. The spectrum of the XYZ Hamiltonian for N = 2m sites in the sector 
with zero momentum coincides with the spectrum of AHq-^ for the twisted staggered 
fermion chain with A^(y) = 3m sites in the sector where T^j^ — (—1)™ after changing 



variables according to ( 39 1 . The multiplicities of the eigenvalues in the two models are 
different. 

4.2 A mapping to hard-particle configurations 

The conjectures [6] and [7] relating the spectra of the XYZ chain and the staggered fermion 
model, raise naturally the question if there is a mapping between the models, at least 
in some subsectors. At the XXZ point C = such a mapping was discussed in p[^: the 
fermion chain with -/V(/) sites and m fermions is equivalent to the twisted spin chain with 
N = N^f-j — m sites and m spins down, the twist being the eigenvalue of the translation 
operator in the fermion model. If we represent an occupied site on the fermion chain by • 
and an empty site by o then the correspondence between fermions and spins is given by 

3 — • — c ^ - and 5 — c ^ + 



This mapping has no direct generalisation to the off-critical case. The reason is the 
absence of conservation of the number of down spins in the general XYZ chain as opposed 
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to the particle conservation in the staggered ferniion chain. However, the path basis was 
designed to implement particle conservation (the number of down steps). Thus, we focus 
on the path states, and try to conceive a mapping between them and the fermion model. 

Let us consider a typical path starting at height £i = i and terminating at some 
^N+i = i + 3p for fixed integer p. As before, let m be the number of decreasing steps. 
Recall that adjacent heights obey £j+i —£j = ±1. However, notice that at rj = 7r/3 we may 
shift any local height variable by a multiple of three without changing the corresponding 



state, as can be seen from the vectors (21 1. Thus, instead of a decreasing step £j+i = — 1 
we can modify the path locally according to £j+i = £j + 2 as shown in figure |5ja). This 
motivates the following construction: given a path we replace each decreasing step by a 
step of two units up, and then continue with usual, appropriately shifted steps up. This 
procedure yields a new, monotone increasing path from height £ to height £ -\- N -\- m -\- 1 as 
illustrated in figure [6|a). Next, we associate to the each of the two types of steps particle 
configurations along the vertical axis according to the rules display in figure [sjb) . This 
is quite reminiscent of the correspondence in the critical case. Thus, we obtain from a 
path a particle configuration with A^(j) — N ra sites and to particles (see figure [6]^a) for 
illustration), with the hard-core rule that particles cannot be adjacent to each other, just 
like in the fermionic case. Notice in particular that because of the condition TV = 2to + 3p, 
the length of the particle chain N(^f^ — 3(to + p) is always a multiple of three. 



^3 



£,+2 / 



(a) (b) 

Figure 5: (a) Local modification of the path, (b) Correspondence between path steps and 
particle configurations. 

For fixed £, the position of the particles • on the vertical axis are given as 

yj^Xj+£ + j, j = l,...,m. 

We fix the origin at ^ = y = 0, and consider periodic boundary conditions, so the j/j's are 
considered mod . 

The mapping between paths and particle configurations is not one-to-one. We illustrate 
this by analysing the effect of a translation on the path configuration by one step to the 
right. The last step of the path is simply removed, and glued to the first one. However, we 
would like to respect the rule that the initial height of the new path is 0, 1 or 2. Thus, a 
vertical shift of all heights by ±3 units might be necessary, and lead to a different particle 
configuration. There are multiple cases, (i) If the last step of the given path goes up, and 
£ — 1,2 the translation has no effect on the configuration of hard-particles, (ii) However, 
if £ = the path has to be shifted by three units, and thus the particle positions are 
cyclically translated by 3: yj — >■ -I- 3 mod TV(/). (Hi) If the last step of the initial path 
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is decreasing, then a translation of the path leaves the particle configuration unchanged if 
1 = 2. (iv) For £ = 0, 1 however, the positions in the particle configuration constructed 
from the translated path are shifted according to yj — >■ ■ + 3 mod Nfjy For given N and 
paths with m steps down, the number of hard-particle configurations obtained through 
the mapping is obtained by counting the paths corresponding to cases (ii) and (iv). This 
yields 



m J \m — 1 J N(f) — TTi\ m 

where we used -/V(j) — N + m. A little combinatorics shows that this is the number 
of possible hard-particle configurations for Af(/) sites and m particles. Also, we see 
that translation of the path configurations is related to translation of the hard-particle 
configurations along the vertical axis by three steps. 

Next, we would like to understand the nature of these particles by examining the local 
action of the XYZ supercharges, explained in section [33] on the particle configurations. 
Recall that the action on the path states corresponds locally to transforming two up steps 
into a single down step. Given our rules identifying paths with particle configurations, it 
is not difficult to see that this corresponds to insertion of a particle while respecting the 
nearest-neighbour exclusion rule (see figure [6]jb)). Recall that this comes with a weight 
(—l)^h(w£^^-^)'^ in terms of the positions for the path. The corresponding position of 
insertion in the hard-particle state is y = x + £ + j where j is the number of particles on 



sites 0, 1, . . . , y — 2. Hence, the weight becomes (— f )^+^ h{wy)'^ as follows from ( 25 1 and 



the periodicity of h{u). The factor (—1)-' in the string suggests that the hard particles 
are indeed fermions. Furthermore it is tempting to use this in order to identify the 
coupling constants Xy. This requires taking into account a systematic identification of the 
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hard-particle states in terms of the path states (what is deHcate as the proposed mapping 
is not one-to-one), changes of normalisations through the supersymmetry operation, and 
finally the restriction to the special momentum spaces for both models. While we are not 
in a position to carry out this program, we nevertheless put forward the following 

Conjecture 8. The coupling constants of the corresponding fermion model are given by 



where Wy is the linear function defined in ( 25 ) 



We see that these coupling constants depend only on the combination s-\-t. The evidence 
for this conjecture is that this parametrisation of the coupling constants uniformises a 
family of elliptic curves appearing in the coordinate direct coordinate Bethe ansatz for the 
fermion chain [53 . In particular, it implies that upon appropriate rescaling the eigenvalues 
of the fermion chain do not depend on s + 



5 Conclusion 

We have studied the XYZ chain and the eight-vertex model along the supersymmetric line, 
and showed that it possesses & M — (2,2) supersymmetry on the lattice. A consequence is 
that chains of different length have common positive energy levels in certain momentum 
sectors, which are organised into supersymmetry quadruplets. Moreover, we presented 
a derivation of the supersymmetry by means of the Bethe ansatz for the eight-vertex 
model, and showed that the supercharges perform simple local operations on the path 
basis. This analysis led us to a novel characterisation for the ground states of the XYZ 
chain with odd length. Finally, we reported some observations that the XYZ chain 
along the supersymmetric line and the staggered supersymmetric fermion chains with 
nearest-neighbour exclusion have exact common eigenvalues in certain subsectors. 

There are many open questions and extensions. To us, it seems most interesting to 
clarify further the nature of the ground states for the chains of odd length. We hope 
that the supersymmetry will be helpful, for instance to prove that there are exactly two 
zero-energy ground states. A central tool in supersymmetric theories is the Witten index 
tr(— 1)^ [28 : it provides a lower bound on the number of zero-energy states. Indeed, 
it would be interesting to define this quantity or find at least a suitable analogue for 
the present theory. As the fermion number coincides with the number of sites, the 
formal generalisation leads to a trace which runs over an infinite collection of Hilbert 
spaces, what one would have to make sense of. Similar considerations apply to the index 
tr((— l)^J^e~'^^) defined in [51]. A possible way to resolve these problems might be to 
establish a more complete mapping between the XYZ chain and the staggered fermion 
chain. For the latter, there exists a standard procedure to find the Witten index, and 
determine the exact number of ground states using cohomology arguments (see e.g. |55jl. 
Further insights into the structural properties of the ground states will certainly be 
obtained by considering the inhomogeneous eight-vertex model, as was the case in the 
trigonometric limit [25] 153] [HI [5^ . Almost all developments in this work considered 
the homogeneous version, and the supersymmetry appears to be intimately related to 
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translation invariance. It would be interesting to see if (and how) this symmetry persists 
in the inhomogeneous case. 

Finally, let us point out that the supersymmetry presented in this article is a particular 
featxire of the 77 = 7r/3 model. It is natural to ask for an extension to general roots-of- unity 
points such as 77 = 7r/(fc + 2) with k = 1,2,3,... The case k = 2 was already addressed in 
|19j from the point of view of fermions with generalised exclusion rules. This allowed the 
identification of a supersymmetric point for the Fateev-Zamolodchikov integrable spin-1 
chain. For more general trigonometric models, the points 77 = TT/{k + 2) were identified as 
the combinatorial points for fused spin-fc/2 models |55], as anticipated in |57j (see also 
Indeed, using the works [SHI EH] we can show that these coincide precisely with the 
cases where a lattice supersymmetry is present. This generalisation will be addressed in a 
forthcoming publication |61j . 

Acknowledgements 

This work was supported by the NSF grants DMR/MPS-0704666 and 1006549. CH was 
supported in part by the ERG AG GONFRA and by the Swiss NSF. He acknowledges 
support from the Kavli Institute for Theoretical Physics under NSF grant PIIY05-51164, 
and would like to thank the KITP, where most of this work was done, for hospitality. 
Moreover, GH would like to thank Luc Blom, Gaetan Borot, Bernard Nienhuis, and 
Paul Zinn- Justin for very fruitful discussions, and Liza Huijse and Anton Zabrodin for 
discussions at early stages of this work. 

A Properties of the supercharges 

In this appendix we present some technical details about the properties of the supercharges 
introduced in section [2] In the first part, we prove the nilpotency property, and in the 
second part, we show that their anticommutator generates the XYZ Hamiltonian. 

A.l Nilpotency 

Let us prove that the operators Qn "have square zero" in the sense that 
Qn+iQn = 0. 

To this end, we need a set of anticommutation rules for the local operators qj defined in 
the main text ([T]). We have the rule 



This can be shown along the lines of [20 , and therefore we only sketch the proof for 
another relation involving q^. Let us first consider q^qj. We find that its action non-zero 
only on states having spins — at position j and N . We find 



qiqj + qj+iqi = 0, I <i < j < N. 



(42) 





• -t- 

N+2 
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Reversing the order of the g's, we have to take into account the shift and therefore consider 
Qj+iQo- Its action yields 

q + + + + +> 

j N ^ j + 1 N j + 1 N 

-C|i--- + +----) + C'l ))• 

j+1 N j + i N ' 

We see that this result coincides with the previous one, except for a minus sign. Combining 
these two equations, we find therefore 

qoqj + qj + iqo = 0, j = 1,...,N -1 (43) 

when acting on 'Hat. 

These relations are useful in order to prove that the supercharges are nilpotent in 



the sense stated above. In a first step, we observe that the (42 1 and (43) can be used to 
reduce the product of the supercharges to 

^ X 1/2 / 



Qn+iQn = ( ( 51 ('^J + l^J + ^j) + lOlN + qN+lQo 



Let us examine the different terms in this sum. The individual terms are non-vanishing 
only if they act on the following states in Hn- 

{q, + iq,+q^)\ ) = C f I ' ' ' - ++ •••)-|--- + +- •••)), J = h---,N 

(<zigo + <Zo')l----) = cfl )-| + ---- + )), (44) 

N \ 1 N+2 1 N+2 J 

(go9Ar + g^v+igo) I • • • -) = (-l)'^+'C f I + • • • - +)-|---- + + 

N \ 1 N+1 1 N+1 

As all the expressions are proportional to C, we see that in the XXZ limit C = the 
relation Qn+iQn = is immediate. Actually, it would not even be necessary to impose 
the restriction to certain momentum spaces in this case. However, for general C the 
relation only survives on the special momentum sectors. Intuitively, this can be seen as 



follows: we see that the operations defined in (44) insert pairs ++ to left and right of a spin 



— , thus we expect that the summation of these on a periodic chain will lead to telescopic 
cancellations. This will however only work in a momentum sector compatible with the 



sign in the third expression in ( 44 ) which coincides with the eigenvalue of the translation 
operator. More concretely, let us consider the case of Qn+iQn on a momentum state 
IV'a) built from spin configuration a = axa^ ■ ■ ■ cxn, that is 

N 
3=0 

We leave aside the issue of normalisation. Let us suppose that a has m spins — and denote 
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their positions by a; i, ... , Xm- Using the rules defined in (44 1 we can write 

N-l 
j=Q 

m 

+ (-1)^+1 (1+) ® T-^' l«) ® 1+) - Tlr'^'la) ® 1++)) 

m 

+ E ® ^^''1") - 1+) ® ^^"'^ l«) ® 1+)) ■ 



where C — \J {N + 2) /N. Because of tff — (—1)^+^ this can be written in terms of the 
positions yi, . . . , UN-m of the individual spins + in a. We find a simplified expression: 



N-7 



CC'QN+iQNli^a) H-^)" E {\+)^T-y^a)(g>\+)^T],~y'\a)(g>\++) 

N-m 

- E (i+- 



>T-y'\a)-\+)^T]fy^\a) 



Now observe that if the configuration a has a spin + at position yi we can write 

\+)^T^y'\a) = Tl,-y^\a)(S\+). 

Using this in the preceding formula we conclude that all terms cancel mutually. This 
proves the statement Qn+iQn = 0. 

A. 2 The Hamiltonian as an anticommutator 

In this appendix we show in detail that if we restrict the Hamiltonian to subsectors where 
the eigenvalue of the translation operator is = (—1)^+^ then it can be written as 
"anticommutator" 

Hn = Q\Qn + Qtv-iQIv-i- (45) 

First, it is useful to introduce the projector on the momentum spaces that we are interested 
in. It is given by 

N-l N-l 

n- = ^E«-^E(-i)'"^^'^n- 

From the definition of the supercharges (|7| and the translation properties established in 
(|8]), we conclude that Qn^n = nAr+iQjv = Qn- Using this, we conclude that 

Q^Qn = N{N +l)nNq\TiN+iqiTiN. 
Qn-iQ^m-i = ^(^ - mNQiIlN-iqlllN- 
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Second, using the definition of the projector, and again the rules (|8|, we find that 

N 



' N-l 



AT- 



We reduce in a third step the sum of these expressions through an application of the 
following anticommutation relations 

q^q]+q]+lq^ = Q, 1 < i < J - 1 < iV - 1, 

qm] + q]+iqQ = o, 2<i<N-i. 



These can be derived in a similar way as the relations ( 42 1 and ( 43 1 . After some algebra, 
we are left with 



njvffArn^l = NUn (q\qi + qiq\ + qUi + qUi) H^^. 
The remaining quadratic terms can be expressed through simple spin operators. We find 

q]g, = i(l + C^)(l-a|), 
(■2 

4. It 

and ^Q^i = q[q2Tj\j. Using these relations and again translation invariance, we conclude 
that 



HmHmU 



NUn (cr^f^s + + C(o-7ct+ + (Ti )) Hjv 



- ivn 



n 



The expression on the right-hand side is nothing but the restriction of the XYZ-Hamiltonian 
(T|2 1 to the momentum sectors with = (—1)^''"^, what proves the statement. 



B Reduction from = 3 to = 2 sites 

In this appendix we show that the supercharge in the path basis can be written as linear 
superposition of the supercharges defined in section |2.2| for the most simple case of three 
and two sites. 

There is a single 7r-momentum state for N = 2 sites. In the canonical basis it is given 
by (we do not normalise the states): 

!</>> = 1+-) -!-+)• 
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For N = 3 sites, there are four possible states which are invariant under translation 



1^3) = 1+—) + + 1— +), iV'4) = I — ). 

Let us now turn to the path basis. For even A'' it is redundant. Indeed, for A'' = 2 

there are i/(2) = 6 admissible paths but the Hilbert space has dimension d = 2^ = 4. 
Indeed, one verifies that the different states are related through the identity 



h{we+i) 



( £ + 1 



\ 



l+l £+1 



( 



h{we-i) 



V 



1 £-1 



e-i I 



In some sense, the relation is trivial here because the difference of the path states on both 
sides is proportional to the singlet state |(/>), and equahty of the proportionality factors on 
both sides is readily verified. 

For iV = 3 there are two states in the path basis which are invariant under translation. 
Up to factors, they are given by 




As mentioned before, the path basis for odd N is incomplete. In this concrete example, 
we see that the two states that are missing have to be invariant under translation. The 
requirement that they are orthogonal to all path states determines them up to linear 
combinations and normalisations. We find it convenient to choose 



IX3)=Cl^l) + l^3>, |X4) = CIV'4) + |V2), (46) 

where we used the coordinate basis. 

In order to find the action on on the vectors we write simply decompose the 
path basis and the two missing states for three sites in terms of the spin basis according 
to IXi) = Ej=i^ijlV'j)- Hence 

4 
i=i 

where the hi are constants (in this example the map is necessarily of rank 1). Hence 

4 

OJlVi) = Y.^A-\ihi\<^). (47) 
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Thus we have to determine the matrix A involved in the change of basis and the vector b. 
For the former it is convenient to abbreviate fj{x) — Y[k=o + 27rfc/3, q^). Then we 
find that 



A- 



( 3h{s) -C/4(s) -C/i(s) 3/4(5) \ 

3/i(i) -Chit) -Chit) 3/4(t) 

C 1 

V 1 c J 



where we used the theta function identity 



{u,q^) (^4 



4 [u,q 



TT 



C^i{ 

-^i(u 



7') ^1 

^4 



(notice that using the definition of C this turns out to be an identity involving products of 
five theta functions and therefore does not simply follow from Riemann's identity). 



Next, let us determine the bj. The only path state which is not annihilated by Q\ is 



Ixi). The application of the local transformation rules defined in (34) and (35 1 leads to 




!+ 1 



-h{we+j 



^2 £ 



Finally the summation over £ = 0, 1, 2 then yields 61 = /i((s - t)/2) J2i=o Kwe)^. From 
the rules for the action of Q2 it is evident that 62 = 0. Not evident however are the values 
of 63 and 64. We follow the proposal made in the main text: the supercharges annihilate 
the two missing states at odd length. Hence we set &3 = 64 = 0. Then it is a simple matter 



to find the action of Q2 on the spin states. After having computed the inverse A of the 



coordinate transformation we find from (47 1 



Ql 



const. X 



{hit)Ql + fiit)Q 



21 . 



with the functions fjit) defined above. The overall factor is a function of s,t and q. The 
action of the operators Q\ and Q\ on the spin states is 



Thus, we see that Q\ — R2Q\Rz where Rn is the spin-reversal operator introduced in 
section [2.1 1 Find with thus the Hermitian conjugates of the supercharges constructed in 
section 
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